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ABSTRACT 


This  project  is  investigating  the  problem  of  locating  seismic  events  from  combined  data  sets  of  regional  and 
teleseismic  arrival  times,  based  on  the  use  of  a  unified  3D  model  of  the  Earth's  velocity  structure  to  predict  travel 
times  for  both  types  of  data.  Inherent  to  this  problem  is  the  joint  tomographic  calibration  of  the  Earth  model  with 
both  regional  and  teleseismic  ground-truth  data.  The  goal  of  the  unified  approach  is  to  remove  the  inconsistencies 
that  result  when  travel-time  predictions  are  performed  with  a  mixture  of  separately  calibrated  regional  and  global 
models  and  empirical  corrections,  which  can  lead  to  degradation  in  event  location  accuracy.  Additionally,  the 
unified  approach  provides  a  rigorous  framework  for  choosing  optimal  relative  weighting  of  different  types  of  data 
used  in  locating  an  event. 

We  are  addressing  a  number  of  practical  difficulties  that  arise  in  the  pursuit  of  a  unified  location/calibration 
capability.  One  is  to  develop  fast  and  accurate  raytracing  techniques  for  modeling  different  types  of  seismic  arrivals. 
A  second  is  the  computational  challenge  of  performing  joint  tomographic  calibration  with  very  large  numbers  of 
data  and  model  parameters.  We  are  considering  various  strategies  for  reducing  the  problem  size,  such  as  averaging 
data  from  proximate  events  and  stations,  and  for  breaking  the  calibration  problem  into  equivalent  sub-problems  by 
region  or  depth  range.  The  calculation  of  model  uncertainty,  and  its  translation  to  travel-time  prediction  uncertainty, 
provides  optimal  data  weighting  for  locating  events,  but  adds  greatly  to  the  computational  challenge. 

A  significant  computational  savings,  in  both  calibration  and  location,  is  achieved  by  employing  a  linearized 
approximation  to  travel-time  calculations,  obtained  by  integrating  the  slowness  function  of  a  3D  model  along  ray 
paths  computed  in  a  ID  reference  model.  While  this  approximation  is  not  adequate  for  regional  travel  times,  we 
have  performed  numerical  experiments  that  indicate  it  is  quite  accurate  for  teleseismic  travel  times,  as  judged  by 
comparing  to  travel  times  calculated  with  a  highly  accurate  raybending  method.  We  have  also  continued  numerical 
experiments,  initiated  last  year,  with  joint  regional/teleseismic  travel-time  tomography  built  on  the  linearized 
forward  model  for  teleseismic  data.  The  results  to  date  indicate  that  biases  between  regional  and  teleseismic 
travel-time  predictions,  which  create  difficulties  for  an  event  location  algorithm,  are  removed  by  the  joint  calibration 
process. 


338 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

SEP  2010 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 


4.  TITLE  AND  SUBTITLE 

A  Unified  Approach  to  Joint  Regional/Teleseismic  Calibration  and  Event 
Location  with  a  3D  Earth  Model 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Massachusetts  Institute  of  Technology, 77  Massachusetts 
Avenue, Cambridge, MA, 02139 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROIECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


13.  SUPPLEMENTARY  NOTES 


Published  in  Proceedings  of  the  2010  Monitoring  Research  Review  -  Ground-Based  Nuclear  Explosion 
Monitoring  Technologies,  21-23  September  2010,  Orlando,  FL.  Volume  I.  Sponsored  by  the  Air  Force 
Research  Laboratory  (AFRL)  and  the  National  Nuclear  Security  Administration  (NNSA).  U.S. 
Government  or  Federal  Rights  License 


14.  ABSTRACT 

This  project  is  investigating  the  problem  of  locating  seismic  events  from  combined  data  sets  of  regional  and 
teleseismic  arrival  times,  based  on  the  use  of  a  unified  3D  model  of  the  Earth’s  velocity  structure  to  predict 
travel  times  for  both  types  of  data.  Inherent  to  this  problem  is  the  joint  tomographic  calibration  of  the 
Earth  model  with  both  regional  and  teleseismic  ground-truth  data.  The  goal  of  the  unified  approach  is  to 
remove  the  inconsistencies  that  result  when  travel-time  predictions  are  performed  with  a  mixture  of 
separately  calibrated  regional  and  global  models  and  empirical  corrections,  which  can  lead  to  degradation 
in  event  location  accuracy.  Additionally,  the  unified  approach  provides  a  rigorous  framework  for  choosing 
optimal  relative  weighting  of  different  types  of  data  used  in  locating  an  event.  We  are  addressing  a  number 
of  practical  difficulties  that  arise  in  the  pursuit  of  a  unified  location/calibration  capability.  One  is  to 
develop  fast  and  accurate  raytracing  techniques  for  modeling  different  types  of  seismic  arrivals.  A  second 
is  the  computational  challenge  of  performing  joint  tomographic  calibration  with  very  large  numbers  of 
data  and  model  parameters.  We  are  considering  various  strategies  for  reducing  the  problem  size,  such  as 
averaging  data  from  proximate  events  and  stations,  and  for  breaking  the  calibration  problem  into 
equivalent  sub-problems  by  region  or  depth  range.  The  calculation  of  model  uncertainty,  and  its 
translation  to  travel-time  prediction  uncertainty  provides  optimal  data  weighting  for  locating  events,  but 
adds  greatly  to  the  computational  challenge.  A  significant  computational  savings,  in  both  calibration  and 
location,  is  achieved  by  employing  a  linearized  approximation  to  travel-time  calculations,  obtained  by 
integrating  the  slowness  function  of  a  3D  model  along  ray  paths  computed  in  a  ID  reference  model.  While 
this  approximation  is  not  adequate  for  regional  travel  times,  we  have  performed  numerical  experiments 
that  indicate  it  is  quite  accurate  for  teleseismic  travel  times,  as  judged  by  comparing  to  travel  times 
calculated  with  a  highly  accurate  raybending  method.  We  have  also  continued  numerical  experiments, 
initiated  last  year,  with  joint  regional/teleseismic  travel-time  tomography  built  on  the  linearized  forward 
model  for  teleseismic  data.  The  results  to  date  indicate  that  biases  between  regional  and  teleseismic 
travel-time  predictions,  which  create  difficulties  for  an  event  location  algorithm,  are  removed  by  the  joint 
calibration  process. 


15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

Same  as 
Report  (SAR) 

18.  NUMBER 
OF  PAGES 

10 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


2010  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


OBJECTIVE 


Previous  studies  have  demonstrated  that  3D  models  of  crust  and  upper  mantle  velocity  structure  predict  significantly 
more  accurate  travel  times  for  regional  seismic  phases  than  do  ID  global  models,  leading  to  more  accurate  locations 
for  small,  regionally  recorded  events  (e.g.,  Ryaboy  et  al.,  2001;  Murphy  et  ah,  2005;  Flanagan  et  ah,  2007;  Reiter 
and  Rodi,  2009).  However,  it  is  not  uncommon  to  also  observe  teleseismic  arrivals  from  small  events  of  monitoring 
interest,  at  least  in  part  of  the  teleseismic  distance  range.  While  adding  teleseismic  to  regional  data  should  in  theory 
improve  event  locations  even  further,  it  has  not  always  been  the  case  in  practice.  For  example,  Yang  et  ah  (2004) 
found  that  event  locations  obtained  by  combining  regional  and  teleseismic  data  were  only  more  accurate  than 
locations  found  from  the  separate  data  sets  if  they  corrected  for  a  bias  between  teleseismic  and  regional  travel-time 
predictions,  each  made  from  a  different  3D  model:  CUB  1.0  (Ritzwoller  et  ah,  2002)  for  regional,  and  J362D28 
(Antolik  et  ah,  2003)  for  teleseismic.  The  correction  they  inferred  is  small  (0.79  s),  but  even  small,  systematic  shifts 
in  teleseismic  travel  times  can  induce  large  location  errors  owing  to  the  high  apparent  velocity  of  teleseismic  waves. 

We  are  investigating  a  methodology  for  seismic  event  location  based  on  the  principle  that  travel  times  for  all  seismic 
phases  should  be  predicted  consistently  with  a  single  3D  Earth  model.  This  unified  modeling  approach  would  be 
applied  as  the  forward  model  in  locating  events  from  combined  sets  of  regional  and  teleseismic  data.  Moreover,  the 
3D  model  would  be  calibrated  with  regional  and  teleseismic  ground-truth  data  in  a  joint  tomographic  analysis.  A 
unified  approach  to  location  and  travel-time  calibration  is  not  a  new  concept  and,  in  fact,  has  been  the  approach 
taken  with  ID  Earth  modeling  for  many  decades.  Applying  the  approach  to  3D  Earth  models,  however,  presents 
many  strategic  and  computational  challenges,  which  this  project  is  attempting  to  address. 

RESEARCH  ACCOMPLISHED 


Our  project  is  focusing  on  the  development  of  a  prototype  capability  for  event  location  and  travel-time  calibration 
with  3D  Earth  models,  applicable  to  regional  and  teleseismic  data.  Our  goal  is  not  to  develop  a  complete  and 
operations-ready  3D  capability  but  rather  to  demonstrate  the  value  and  practicality  of  the  unified  approach.  The 
following  sections  describe  our  progress  to  date  in  developing  two  key  components  of  the  approach:  3D  travel-time 
calculation  and  joint  regional/teleseismic  tomography. 

Travel-Time  Linearization 

A  number  of  numerical  methods  are  available  for  raytracing  and  travel-time  calculation  in  3D  Earth  models,  such  as 
the  finite-difference  eikonal  method  (e.g.,  Podvin  and  Lecomte,  1991),  fast  marching  (Rawlinson  and  Sambridge, 
2004),  and  two-point  raybending  (e.g.,  Um  and  Thurber,  1987).  Given  the  computational  intensity  of  these  hill  3D 
methods,  many  applications,  such  as  global  tomography  with  teleseismic  data,  have  used  approximate  travel  times 
calculated  with  linearization  around  a  ID  reference  model.  In  this  approximation,  the  travel  time  T  through  a  3D 
model  of  interest  is  obtained  by  integrating  the  model's  slowness  function  along  the  Fermat  raypath  computed  in  a 
different,  reference  model  velocity  model: 

T  «  /  dx  a(x;  izref)  v'1  (x)  (1) 

where  x  is  the  position  vector;  v(x)  is  the  3D  velocity  function  of  interest;  and  a(x;  vref)  is  a  sensitivity  distribution 
concentrated  on  the  Fermat  raypath  evaluated  in  the  reference  velocity  function,  vret{x).  When  vref  =  v,  T  is  exact. 
When  vref  is  ID,  linearization  of  travel  times  offers  great  computational  advantages  in  both  computation  time  and 
storage.  It  also  greatly  simplifies  the  computation  of  times  for  secondary  arrivals. 

Rodi  et  al.  (2009)  conducted  a  preliminary  investigation  of  the  accuracy  of  travel-time  linearization  in  3D  velocity 
models.  Their  experiments  involved  the  computation  of  P-wave  travel  times  for  various  paths  in  southern  Asia 
through  the  JWM  model  of  Reiter  and  Rodi  (2009).  Linearized  times,  using  reference  rays  for  model  AK135 
(Kennett  et  al.,  1995),  were  compared  to  travel  times  calculated  with  the  Podvin-Lecomte  (P-L)  finite-difference 
method,  with  the  presumption  that  the  latter  are  more  accurate  and  can  be  used  to  judge  the  adequacy  of 
linearization.  The  numerical  tests  reported  last  year  showed  that  linearization  about  AK135  was  very  inaccurate  for 
computing  travel  times  at  regional  distances,  as  expected,  but  seemed  to  be  a  good  approximation  at  teleseismic 
distances  except  for  some  bothersome  artifacts  related  to  shallow  velocity  structure.  The  tests  were  not  definitive, 
however,  because  the  accuracy  of  the 
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P-L  times,  which  are  also  approximate,  was  not  known.  Moreover,  the  JWM  model  is  the  same  as  AK135  below  a 
depth  of  410  km,  weakening  the  tests. 

We  have  continued  this  investigation  of  travel-time  linearization  with  the  aim  of  getting  a  more  definitive 
assessment  of  its  accuracy.  The  new  tests  we  report  here  involve  three  main  changes  from  our  preliminary  tests. 
First,  the  calculations  were  performed  with  the  3D  velocity  model  taken  to  be  a  portion  of  Lawrence  Livermore 
National  Laboratory’s  (LLNL's)  recently  developed  G3D-LLNL  global  tomographic  model  (Simmons  et  al.,  2010), 
which  has  3D  velocity  anomalies  throughout  the  crust  and  mantle.  Second,  we  modified  our  algorithm  for 
computing  linearized  travel  times  to  be  more  accurate  when  reference  rays  traverse  low  velocity,  shallow  layers  in 
the  3D  model. 

In  a  third  departure  from  our  earlier  investigation,  we  used  highly  accurate  raybending  calculations  as  the 
benchmark  for  assessing  the  accuracy  of  linearization.  The  raybending  calculations  through  G3D-LLNL  were 
performed  by  LLNL  using  an  algorithm  they  adapted  from  the  3D  raytracing  procedure  of  Zhao  (1992)  and  Zhao 
and  Lei  (2004).  This  technique  combines  the  pseudo-bending  method  for  continuous  media,  developed  by  Um  and 
Thurber  (1987),  and  refraction  (Snell's  law)  at  velocity  discontinuities.  The  specific  approach  is  outlined  in 
Simmons  et  al.  (2009). 

Our  new  tests  considered  travel  times  between  various  station  locations  and  a  3D  grid  of  event  locations.  The  event 
grid  covers  the  box  15°-55°  N,  30°-80°  E  (roughly,  the  Middle  East)  at  1°  spacing,  but  contains  only  a  few  event 
depths  between  the  surface  and  120  km.  We  performed  calculations  for  four  station  locations,  one  at  each  corner  of 
the  event  grid  and  each  at  zero  depth.  We  will  refer  to  these  as  Station  NW,  etc. 

Figure  1  compares  travel  times  calculated  through  G3D-LLNL  (with  raybending)  to  those  calculated  through 
AK135.  We  see  that  the  effects  of  3D  velocity  structure  in  G3D-LLNL  on  regional  times  is  several  seconds.  The 
3D  effects  on  teleseismic  travel  times  (distance  greater  than  ~23  deg)  are  smaller  but  still  significant.  For  the  same 
stations  and  grid  of  surface  events,  Figure  2  shows  the  discrepancy  between  travel  times  calculated  through  G3D- 
LLNL  with  raybending  and  with  linearization  around  AK135,  i.e.  integrating  the  slowness  function  of  G3D-LLNL 
along  rays  generated  in  AK135.  Given  the  high  accuracy  of  the  raybending  method,  the  travel-time  discrepancy 
plotted  in  Figure  2  can  be  interpreted  as  the  error  in  the  linearization  calculations.  The  same  results  are  shown 
geographically  in  Figure  3. 

Figures  2  and  3  demonstrate  that  the  errors  in  linearized  Pn  times  can  be  very  large,  as  much  as  30  s  for  some  paths 
(not  seen  in  the  plots  due  to  clipping).  The  errors  at  teleseismic  distances  are  much  smaller,  being  generally  well 
below  0.5  s.  Flowever,  it  is  notable  in  Figure  3  that  the  largest  errors  in  linearized  times  for  these  surface  events 
correlate  with  the  areas  having  the  lowest  near-surface  velocities,  i.e.  the  oceans.  For  deeper  events,  the  surface- 
related  artifacts  disappear  (Figure  4).  While  in  practice  we  do  not  need  to  be  concerned  with  event  locations  in 
oceans  and  slow  sediments,  we  are  still  investigating  whether  our  ray  integration  scheme  may  be  undersampling 
thin,  low  velocity  layers  in  the  3D  model,  or  whether  we  have  simply  reached  a  situation  in  which  linearization  is 
not  adequate. 

Figures  2-4  also  show  that  the  errors  in  linearized  teleseismic  times  are  predominantly  negative,  which  seemingly 
violates  Fermat's  Principle.  One  explanation  for  the  negative  bias  is  a  quirk  in  our  experiments  whereby  the  effects 
of  the  Earth's  ellipticity  were  calculated  differently  in  the  raybending  and  linearization  methods.  We  may  also  be 
seeing  the  above-mentioned  problems  with  thin,  low-velocity  layers  beneath  the  stations  (rather  than  the  events). 
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Figure  1.  Difference  between  travel  times  through  models  G3D-LLNL  and  AK135,  calculated  between  each  of 
four  stations  and  a  grid  of  surface  events.  The  G3D-LLNL  times  were  calculated  with  raybending; 
AK135  times  were  calculated  analytically.  The  time  differences  are  plotted  as  a  function  of 
epicentral  distance. 


Station  NW  Station  NE 


Figure  2.  Discrepancy  between  travel  times  calculated  through  model  G3D-FFNF  with  two  methods: 

raybending  and  linearization  relative  to  AK135  rays.  Each  panel  shows  the  discrepancy  vs.  distance 
for  paths  from  one  of  the  corner  stations  to  the  geographic  grid  of  surface  events.  The  travel-time 
scale  has  been  clipped  at  8  sec. 
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Joint  Regional/Teleseismic  Tomography 

Rodi  et  al.  (2009)  reported  on  preliminary  tests  of  joint  regional/teleseismic  calibration,  whereby  the  JWM  model  of 
southern  Asia  derived  by  Reiter  and  Rodi  (2009)  was  updated  through  linearized  tomography  applied  to  a 
combination  of  regional  and  teleseismic  P-wave  travel  times.  Linearized  tomography  uses  Equation  (1)  to  model 
travel  times.  For  teleseismic  times,  we  take  the  reference  velocity  model  to  be  AK135,  which  our  numerical 
experiments  have  shown  to  be  adequate.  The  same  experiments  showed  that  linearization  about  AK135  is  not 
adequate  for  regional  travel-time  calculation;  therefore,  we  used  JWM  as  the  reference  model  for  regional  times. 

The  regional  data  set  used  in  the  joint  calibration  tests  was  identical  to  that  used  by  Reiter  and  Rodi  (2009)  in 
constructing  JWM.  Two  teleseismic  data  sets  were  considered,  both  extracted  from  the  EFIB  database  (Engdahl  et 
al.,  1998)  and  both  including  arrivals  only  in  the  epicentral  distance  range  35°-65°.  In  the  first  teleseismic  data  set 
(termed  TelE),  the  stations  were  restricted  to  be  located  in  or  near  the  JWM  study  region,  i.e.,  in  the  box  0°-60°  N, 
30°-120°  E.  The  event  locations  were  unrestricted  (except  for  the  distance  criterion)  and  fell  primarily  outside  this 
box.  The  second  data  set  (TelS)  restricted  the  events  to  be  within  the  study  box  and  accepted  all  stations  obeying  the 
distance  criterion.  The  regional  and  both  teleseismic  data  sets  included  only  events  which  have  EHB  focal  depths 
less  than  200  km  and  which  are  deemed  well-located  in  accordance  with  criteria  described  by  Reiter  and  Rodi 
(2009).  To  reduce  the  size  and  redundancy  of  the  data  sets,  individual  events  were  lumped  into  summary  events  in 
the  same  manner  as  done  for  JWM,  with  each  summary  event  being  a  node  on  a  0.5°  geographic  grid  and  an 
irregular  depth  grid. 


Lin.  minus  R-B  Time  (s) 


-1  -0.5  0  05  1 


Figure  3.  Discrepancy  between  travel  times  calculated  through  model  G3D-LLNL  with  raybending  and 

linearization.  The  discrepancies  shown  are  the  same  as  in  Figure  2  but  plotted  as  a  function  of  the 
epicenter  of  the  surface  events.  The  travel-time  scale  has  been  clipped  at  -1  and  +1  sec. 
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We  have  repeated  these  experiments  using  our  improved  algorithm  for  computing  linearized  travel  times,  and  with 
upgraded  teleseismic  data  sets.  The  teleseismic  data  were  extracted  from  the  updated  EHB  database  distributed  by 
the  International  Seismological  Centre  (http://www.isc.ac.uk/EHB/index.html).  The  new  database  included  two 
additional  years  (2005  and  2006)  as  well  as  additional  arrivals  for  years  prior  to  2005.  After  forming  summary 
events,  data  set  TelE  contained  41 1,232  teleseismic  P  arrivals  for  15,576  summary  events  and  900  stations.  TelS 
contained  360,706  arrivals  from  4,253  summary  events  and  2,846  stations.  The  teleseismic  data  sets  are 
considerably  larger  than  the  regional  data  set,  which  contains  104,065  arrivals  (>99%  being  the  Pn  phase)  from 
3,689  summary  events  and  603  stations.  These  data-set  sizes  approach  the  capacity  of  our  current  tomography 
software,  which  is  why  we  considered  teleseismic  events  and  stations  separately,  and  why  we  capped  epicentral 
distances  at  65°. 
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Figure  4.  Discrepancy  between  travel  times  calculated  through  model  G3D-LLNL  with  raybending  and 

linearization.  Each  panel  shows  the  discrepancy  for  paths  from  Station  NE  to  a  geographic  grid  of 
events  at  the  labeled  focal  depth  (z). 

Figure  5  shows  a  vertical  slice  through  each  of  two  joint  regional/teleseismic  tomography  models,  corresponding  to 
the  two  versions  of  the  teleseismic  data  set  that  were  combined  with  the  common  regional  data  set.  The  model 
obtained  with  TelE  as  the  teleseismic  data  set  is  labeled  JWMTelE,  and  the  other  model  as  JWMTelS.  Both 
models  are  displayed  as  deviations  from  JWM.  We  see  that  the  models  are  quite  similar  to  each  other,  and  overall 
have  lower  velocities  than  JWM.  The  largest  differences  from  JWM  occur  at  depths  below  200  km,  where  the 
P  velocity  is  significantly  slower. 

The  reason  for  the  lower  mantle  velocities  in  JWM  TelE  and  JWM  TelS  is  clear  if  we  compare  their  travel-time 
residuals  to  those  of  JWM.  Figure  6  displays  separate  histograms  of  regional  (left)  and  teleseismic  (right)  residuals 
for  each  joint  inversion  model  (red  lines),  compared  in  each  case  to  the  corresponding  histograms  for  JWM 
(blue  lines).  We  see  that  JWM  yields  regional  residuals  with  a  mean  near  zero  but  teleseismic  residuals  with  a  large 
positive  mean:  1.8  s  for  the  TelE  residuals  and  1.6  s  for  the  TelS  residuals.  Both  joint  inversion  models  maintain  the 
unbiased  fit  to  the  regional  data  and  yields  zero-mean  teleseismic  residuals  as  well.  Each  joint  regional  tomography 
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removed  the  teleseismic  travel-time  bias  of  JWM  by  lowering  the  P  velocity  below  200  km  depth,  in  a  way  that 
preserved  the  fit  to  the  Pn  data. 
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Figure  5.  A  vertical  section  through  two  3D  P  velocity  models  obtained  by  joint  inversion  of  regional  and 

teleseismic  P-wave  travel  times,  each  shown  in  terms  of  its  velocity  deviation  from  model  JWM.  The 
top  model  (JWMTelE)  used  teleseismic  data  from  local  stations  and  teleseismic  events,  while  the 
bottom  model  (JWMTelS)  used  teleseismic  data  from  local  events  and  teleseismic  stations.  The 
regional  data  were  the  same  in  both  cases  and  are  from  stations  and  events  inside  the  model  region. 
Note  that  the  color  scale  for  velocity  deviation  is  asymmetrical  (shifted  toward  negative  deviations). 
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Figure  6.  Red  lines:  Distributions  of  regional  (left)  and  teleseismic  (right)  travel-time  residuals  for  models 
JWM  TelE  (top  panels)  and  model  JWM  TelS  (bottom  panels).  The  teleseismic  data  set  is  TelE 
(teleseismic  events)  for  the  top-right  panel  and  TelS  for  the  bottom-right  panel.  Blue  lines: 
Analogous  distributions  for  model  JWM. 
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Joint  Tomography  with  Hypocenter  Relocation 

The  above  tomography  results  and  those  presented  in  Rodi  et  al.  (2009)  do  not  account  for  the  effect  of  event 
mislocations  on  the  velocity  model.  We  have  addressed  this  by  extending  our  tomography  algorithm  to  allow  event 
hypocenters  and  origin  times  to  vary,  in  accordance  with  specified  prior  errors.  We  followed  the  approach  developed 
by  Spencer  and  Gubbins  (1980)  and  Pavlis  and  Booker  (1980),  whereby  a  projection  operator  is  applied  to  the 
travel-time  residual  vector  and  its  sensitivity  matrix  to  remove  patterns  in  the  residuals  that  can  be  explained  by 
perturbing  the  hypocenters  and  origin  times  of  the  events.  Departing  from  the  original  formulations  of  the  approach 
in  terms  of  matrix  inversion,  we  implemented  the  projection  operations  within  the  conjugate  gradients  technique  we 
use  for  computing  tomographic  models. 


To  test  the  new  capability,  we  repeated  the  two  joint  regional/teleseismic  tomography  experiments  of  the  previous 
section,  allowing  variable  event  locations  with  a  prior  hypocentral  uncertainty  of  10  km  (epicenter  and  depth)  and 
origin  time  uncertainty  of  2  s  (all  at  95%  confidence).  Figure  7  shows  a  slice  through  the  resulting  inversion  models 
(called  JWMTelEloc  and  JWMTelEloc).  We  see  that  model  JWMTelEloc  deviates  much  less  from  JWM 
than  did  model  JWMTelE  (Figure  5  top  vs.  Figure  7  top),  showing  that  perturbations  to  the  teleseismic  event 
locations  can  fit  a  significant  portion  of  the  JWM  residuals.  Models  JWMTelSloc  and  JWMTelS  are  less 
different.  However,  as  before,  both  joint  inversion  models  reconcile  the  regional/teleseismic  travel-time  bias  of 
JWM,  as  seen  in  Figure  8.  We  also  see  that  the  residual  histograms  for  JWM  TelE  loc,  and  to  a  lesser  extent 
JWM  TelS  loc,  are  more  peaked  than  before  as  a  result  of  allowing  perturbations  in  the  event  locations. 
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Figure  7.  A  vertical  section  through  two  3D  P  velocity  models  (JWM  TelE  loc  and  JWM  TelS  loc)  obtained 
by  joint  inversion  of  regional  and  teleseismic  P-wave  travel  times,  each  shown  in  terms  of  its  velocity 
deviation  from  model  JWM.  This  figure  is  analogous  to  Figure  5  except  that  the  joint  inversion 
models  were  obtained  in  conjunction  with  linearized  hypocenter  relocation. 
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Figure  8.  Red  lines:  Distributions  of  regional  (left)  and  teleseismic  (right)  travel-time  residuals  for  models 
JWM  TelE  loc  (top  panels)  and  model  JWM  TelS  loc  (bottom  panels).  Blue  lines:  Analogous 
distributions  for  model  JWM. 


CONCLUSIONS  AND  RECOMMENDATIONS 


We  continue  to  develop  and  improve  the  tools  needed  to  perform  event  location  and  travel-time  calibration  with 
regional  and  teleseismic  arrival-time  data,  based  on  a  unified  3D  Earth  model  for  travel-time  calculation.  Our  latest 
numerical  raytracing  experiments  show  that  the  errors  in  linearized  teleseismic  travel-time  calculations  are  0.5  s  or 
less,  which  is  adequate  for  location  and  calibration  but  large  enough  to  warrant  further  effort  on  making  the  errors 
even  smaller.  The  use  of  linearization  for  travel-time  prediction  vastly  reduces  the  computational  requirements  of  a 
unified  approach. 

Our  latest  tests  of  joint  regional/teleseismic  tomography  confirm  the  value  of  this  approach  in  generating  a  3D  Earth 
model  that  predicts  consistent  travel  times  for  regional  and  teleseismic  P  waves.  The  tests  also  show  that  accounting 
for  event  mislocations  in  the  tomographic  procedure  does  affect  the  velocity  model  that  results,  implying  that 
location  uncertainty  must  be  quantified  accurately  in  order  to  properly  resolve  the  location/velocity  trade-off.  The 
use  of  secondary  seismic  phases  (e.g.  surface  reflections  and  S  phases)  may  ultimately  address  this  problem.  We  are 
attempting  to  increase  the  capacity  of  our  tomography  algorithm  significantly  in  order  to  handle  much  larger, 
multiple-phase  data  sets  in  a  desktop  computing  environment. 
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